Nonautonomous discrete rogue waves and interaction in the generalized 
Ablowitz-Ladik-Hirota lattice with variable coefficients 
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We analytically investigate the nonautonomous discrete rogue wave solutions and their interaction 
in the generalized Ablowitz-Ladik-Hirota lattice with variable coefficients, which possess complicated 
wave propagations in time and are beyond the usual discrete rogue waves. When the amplitude of 
the tunnel coupling coefficient between sites decreases, these nonautonomous discrete rogue wave 
solutions become localized in time after they propagate over some certain large critical values. 
Moreover, we find that the interaction between nonautonomous discrete rogue waves is elastic. In 
particular, these results can reduce to the usual discrete rogue wave solutions when the gain or loss 
term is ignored. 
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I. INTRODUCTION 

Rogue waves (alias as freak waves, monster waves, 
killer waves, giant waves or extreme waves), as an im- 
portant physical phenomenon, are localized both in space 
and time and depict a unique event that 'appears from 
nowhere and disappears without a trace' Rogue 
waves (RWs) are also known as 'rogons' if they reap- 
pear virtually unaffected in size or shape shortly after 
their interactions 0]. The study of RWs has become a 
significant subject in many fields since they can signal 
catastrophic phenomena such as thunderstorms, earth- 
quakes, and hurricanes. RWs have been found in the 
ocean [a, 0], nonlinear optics 0-0], Bose-Einstein conden- 
sates [1], the atmosphere 0, and even the finance [Tol | . 
Moreover, some experimental observations have shown 
that optical RWs do exist and play a positive role in non- 
linear optical fibres [B-0|- Such optical RWs differ from 
oceanic RWs that play a negative role and lead to many 
accounts of such waves hitting passenger ships, container 
ships, oil tankers, fishing boats, and offshore and coastal 
structures, sometimes with catastrophic consequences [3] . 
In particular, the analytical RWs have been obtained for 
the nonlinear Schrodinger (NLS) equation [l|, [ll|, G3 , as 
well as some of their extensions with the varying coef- 
ficient s B , the higher orders or the higher dimen- 
sions [14j . 

The discrete NLS equation [l5| and the Ablowitz-Ladik 
(AL) lattice 0, H3|, as two prototypical discretizations 
of the continuum NLS equation, have been studied ex- 
tensively in the field of nonlinear science. The former 
is nonintegrable, but has some interesting applications 
of physics fla - l^ . The latter is integrable an d p os- 
sesses an infinite number of conservation laws [l6l.ll7|. as 
well as has been depicted as an effective lattice to study 
properties of the intrinsic localized modes (23[. More- 
over the nonintegrable discrete NLS equation can also 



be regarded as a perturbation of the integrable AL lat- 
tice [2J]. In addition, the Salerno model (SM) has also 
been presented, interpolating between the nonintegrable 
discrete NLS equation (/x = 0) and the intregable AL 
lattice (e = 0) in the form [25l - t30| 

#„,t + (l/'n+l-|-'0n-l)(l+MlV'n| 2 )-|- £ | | 2 VVi + V1p n = 0, (1) 

which can be derived on the basis of a variational prin- 
ciple SCsM/^n = from the Lagrangian density 

n 

tp n -lW n \^n\ 2 + e|^ n | 4 + 2*#n| 2 , (2) 

where ip n = ip n (t) stands for the complex field ampli- 
tude at the nth site of the lattice, the parameter /i is 
the intersite nonlinearity and corresponds to the nonlin- 
ear coupling between nearest neighbors, e measures the 
intrinsic onsite nonlinearity, and v describes the inhomo- 
geneous frequency shift. The SM has been applied in bi- 
ology [25[ , and Bose-Einstein condensates [3(| • Recently, 
the discrete RWs have also drawn much attention. The 
discrete NLS equation (i.e., for the case fi = 0, e = 1, 
and v = — 2 in Eq. ([T]))has numerically been verified to 
support discrete RWs [31(. The SM has also been found 
to admit discrete RWs from the viewpoint of statistical 
analysis [32|. More recently, it has also been shown that 
exact discrete RWs [H, [34| can exist in the AL lattice 
(i.e., for the case fx = 1, e = 0, and v = —2 in Eq. (Q])) 
on the basis of the limit cases of their multi-soliton so- 
lutions [351 ] . However, there is no reports to date about 
nonautonomous discrete RWs except for the AL lattice 
and the discrete Hirota equation [H, [34| . 

In the present paper, we will explore exact nonau- 
tonomous discrete RW solutions and interaction of the 
generalized Ablowitz-Ladik-Hirota (ALH) lattice with 
variable coefficients given by Eq. ([3]), i.e., the generalized 
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case of Eq. ([T]) without the intrinsic onsite nonlinear- 
ity, where the tunnel coupling coefficients and the inter- 
site nonlinearity are the time-modulated, complex- valued 
and real-valued functions, respectively, the inhomoge- 
neous frequency shift are space- and time-modulated, 
real-valued functions, and the time-dependent gain or 
loss term is added. To do so, we will make use of the 
differential-difference symmetry analysis that can con- 
nect this equation with variable coefficients with the sim- 
pler ones. We show that, for the attractive intersite non- 
linearity, the generalized ALH lattice with variable co- 
efficients can support nonautonomous discrete RWs in 
terms of the rogue wave solutions of the discrete Hirota 
equation. Moreover, we also exhibit wave propagations 
of nonautonomous discrete RW solutions and their inter- 
action for some chosen parameters and functions. 

The rest of this paper is organized as follows. In Sec. 
II, we introduce the generalized ALH lattice with variable 
coefficients, which contains some special lattice models, 
such as the AL lattice, the discrete Hirota equation, and 
the generalized AL lattice. In Sec. Ill, by analyzing the 
phase of the complex field amplitude, we systematically 
present a similarity transformation reducing the general- 
ized ALH lattice with variable coefficients to the discrete 
Hirota equation. In Sec. IV, we determine the self-similar 
variables and constraints satisfied by the coefficients in 
Eq. Moreover, we analyze the relations among these 
functions such that we find that the intersite nonlinear- 
ity (the inhomogeneous frequency shift) is related to the 
gain or loss term (the tunnel coupling coefficient). Sec. V 
mainly discusses nonautonomous discrete RW solutions 
and their interaction of Eq. ((3]) for some chosen param- 
eters and functions. For the given periodic gain or loss 
term, when the amplitude of the tunnel coupling coeffi- 
cient between sites decreases, these nonautonomous dis- 
crete rogue wave solutions are localized in space and keep 
the localization longer in time, which differ from the usual 
discrete rogue waves. Finally, we give some conclusions 
in Sec. VI. 



II. THE GENERALIZED 
ABLOWITZ-LADIK-HIROTA LATTICE WITH 
VARIABLE COEFFICIENTS 

We here address the generalized Ablowitz-Ladik- 
Hirota (ALH) lattice with variable coefficients modeled 
by the following lattice 

M ntt + [A(t)*„+i + A*(i)*„-i] [1 + s(t)|*„| 2 ] 

(3) 

-2«„(t)#„ + i 7 (i)#„ = 0, 

which can be derived in terms of a variational principle 
5Calh/S^'* 1 = from the following Lagrangian density 

£alh = E„ i(K*n,t - * n n,t) + 4 Rc [A(t)*;*„+i] 

+ 5 (t)[A(f)*„+i+A*(t)* n _ 1 ]tf*|* n | a (4) 

-2[2«„(t)-; 7 (t)]|*„| 2 , 



where ty n = ^ n (t) stands for the complex field amplitude 
at the nth site of the lattice, the complex-valued func- 
tion A(t) is the coefficient of tunnel coupling between 
sites and can be rewritten as A(t) = a(t) + ij3(t) with 
a(t) and j3(t) being differentiable, real-valued functions, 
git) stands for the time-modulated intersite nonlinearity, 
v n (t) is the space- and time- modulated inhomogeneous 
frequency shift, and "fit) denotes the time-modulated ef- 
fective gain or loss term. 

In fact, this nonlinear lattice model ([3]) contains many 
special lattice models, such as the AL lattice for the 
case a(t) = const.. Pi t) = v„(t) = y(t) = 0, and 
git) = const. 0, [lT]], the AL equation with addi- 
tional term accounding for dissaption for the case a(t) = 
const., Pit ) = v n (t) = 0, j(t) — const., and g(t) = 
const. [26[, the discrete Hirota equation for the case 
a(t) = const.. Pit) — const., v n (t) = 7(f) = 0, and 
g(t) = const. [3a], the generalized AL lattice given by 
Eq. (TT]) for the case a(t) = const., P(t) = j(t) = 0, and 
git) = const. and the discrete modified KdV equa- 
tion for the case ait) = v„(t) — 7(f) — 0, P(t) — const., 
and g(t) = const. [37| . 

III. DIFFERENTIAL-DIFFERENCE 
SIMILARITY REDUCTIONS AND 
CONSTRAINT EQUATIONS 

We consider the spatially localized solutions of Eq. , 
i.e., limi^i^Qo j^^i)! = 0. To this aim, we search for a 
proper similarity transformation connecting solutions of 
Eq. ([3]) with those of the follo wing discrete Hirota equa- 
tion with constant coefficients [35| , namely 

l $„, r +(A$„ +1 + A*$„_ 1 )(l + |$„| 2 )- 2Re ( A ) $ «= ' ( 5 ) 
which can be derived in terms of a variational principle 
SCii/S^n = from the following Lagrangian density 

Aj = E„ i(^$«,r-*n*Ji, T ) + 4 Re (A$;$ n+1 ) 

(6) 

+ (A$ n+1 +A*$„_ 1 )$ I *J$„| 2 -4Rc(A)|$ n | 2 , 

where $ n = ^ n (r) is a complex dynamical variable at 
the nth site of the lattice, r = r(t) is a real- valued func- 
tion of time to be determined, and the complex-valued 
parameter A can be rewritten as A = a + ib with a and b 
being real- valued parameters. The discrete Hirota model 
§5§ contains some special physical models, such as the 
AL lattice for the case a = 1 and 6 = [lH, [13] and 
the discrete mKdV equation for the case a = and 
6=1 37| ■ It has been shown in Ref. [38[ that the discrete 
Hirota equation ([5|) is in fact an integrable discretization 
of the three-order NLS equation (also known as the Hi- 
rota equation) [39| 

iq t + a(q xx + \q\ 2 q) - ib{q xxx + 6 1 ^ | ^ ) = 0, (7) 

which plays an important role in nonlinear optics |40ll4l|. 

To show the above-mentioned aim, we here apply the 
similarity transformation in the form 

tt n (f) = p(t)e^-W*„[r(*)] (8) 
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to Eq. ([3]), where the function pit) and phase (p n (t) are 
both real-valued functions of indicated variables to be 
determined. To conveniently substitute ansatz ([8]) into 
Eq. and to further balance the phases in every term 
in Eq. ([3]), i.e., ^ n+1 (t), and *„_i(t), we should 

firstly know the explicit expression of the phase tp n (t) in 
transformation ([5]) in space. 

Here wc consider the case that the phase is expressed 
as a quadratic polynomial in space with coefficients being 
functions of time in the form (p n (t) = P2(t)n 2 + p\{t)n + 
Po(t) with po,i,2(i) being functions of time, which is sim- 
ilar to the phases in the continuous NLS (or GP) equa- 
tions with variable coefficients [42j ]. Based on the sym- 
metry analysis, we balance the coefficients of these terms 
^n+i{t), ^n(t), and ty n -i(t) such that we find that the 
phase in transformation ([SJ should be a first degree poly- 
nomial in space with coefficients being functions of time, 
namely 



fn(t) = pi(t)n+p (t), 



(9) 



where po t i(t) are functions of time to be determined. 

Eq. ([5]) with the condition © allows us to reduce 
Eq. ([3]) to Eq. (|5|), variables in this reduction can be de- 
termined from the requirement for the new complex field 
amplitude $ n (r(t)) to satisfy Eq. Thus, we substi- 
tute transformation ([5]) into Eq. ^ along with Eq. © 
and after relatively simple algebra obtain the following 
system of ordinary differential equations 

/>(*) + 7(t)p(*)=0, (10a) 
af(t) - a(t) cosjJx(i) + /3(t) sm Pl (t) = 0, (10b) 
[bP(t) + aa(t)} sinpi(t) + [a/3(t) - ba(t)} 

x cospi(t) = 0, (10c) 
2v n (t) + pi(t)n+p (t) + 2[/3(t) sinpi(t) 

—a(t) cospi(i)] = 0, 



g(t)p 2 (t) = 1, 



(lOd) 
(10e) 



where the dot denotes the derivative with respect to time. 

Therefore, if system (fpQ|) is consistent, then we have 
constructed an algorithm generating nonautonomous so- 
lutions of Eq. (j3)) based on transformation (|8|) and solu- 
tions of Eq. (0: 

Firstly, we solve Eqs. (|10a[) - (|10c[) to obtain the func- 
tions p(t), r(t), and pi(t) in transformation (|8|). 

And then wc consider Eqs. (|10d|) and (|10c[) to deter- 
mine the inhomogencous frequency shift v n (t) and the 
intcrsite nonlincarity g(t) in Eq. ((3j> in terms of above- 
obtained functions pit), r(<), andpi(i). 

Thus, wc have established a similarity transforma- 
tion ([5]) connecting solutions of Eq. ([5]) and those of 
Eq. ([3]). In particular, we here exhibit our approach in 
terms of two lowest-order discrete rogue wave solutions 
of Eq. ([5]) as seeding solutions to find nonautonomous 
discrete rogue wave solutions of Eq. © . 



IV. DETERMINING THE SIMILARITY 
TRANSFORMATION AND COEFFICIENTS 



It follows from Eqs. (|10a|) - (|10c|) that we can obtain the 
variables pit), pi{t) and r(t) in transformation (|8]) in the 
form 



p(t) = p exp 
P\(t) = tan -1 



7(s)ds 



ba(t) - a/3(t) 
aa(t) + bp(t) 



(11a) 
(lib) 



r(t) = (a 2 +b 2 )- 1/2 / [ a 2 {s)+P 2 { S )] 1,2 ds, (11c) 
Jo 

where po is an integration constant. Now it follows from 
Eqs. (|10d|) and (flOe)) along with Eqs. (fTTa|) and (|llb[) 
that we further find the inhomogencous frequency shift 
v n (t) and intersite nonlincarity g(t) in the form 



9{t) = Po 2 exp 



7(s) ds 



v„(t) = vi(t)n + v (t), 



(12a) 
(12b) 



where we have introduced two functions in the inhomo- 
geneous frequency shift v n (t) in the form 



«i(t) = 



vo(t) 



a(t)$(t) - a(t)P(t) 
2[a 2 (t) + p{t)\ ' 

a 2 it)+/3 2 (ty 1/2 



b 2 



Po{t) 



(13a) 
(13b) 



where po(t) is an arbitrary diffcrcntiablc function of time. 

It follows from Eqs. (|12a|) - (|13b|) that, in these coeffi- 
cients of Eq. ([3]), the intersite nonlincarity g(t) (the inho- 
mogeneous frequency shift v n (t)) is related to the gain or 
loss term 7(4) (the tunnel coupling A(i) ~ a(t) + i/3(i)). 
This means that only two varying coefficients (e.g., j(t) 
and A(i) = a(t)+i/3(t)) are left free. Moreover, it follows 
from Eq. ()12a|) that the intersite nonlincarity g(t) is al- 
ways positive (i.e., the attractive intersite nonlincarity). 
In addition, it follows from Eq. (jllaj) that the gain or 
loss term 7(f) can also control the function pit), which 
is used to modulate the amplitude of the complex field 

#»(*)■ 

For the inhomogeneous frequency shift w„(i) given by 
Eq. (|12bj) . when a(t) 7^ cj3{t) with c being a constant, 
the inhomogeneous frequency shift v n (t) is a linear func- 
tion of space n with coefficients being functions of time. 
In the absence of the discrete space n in the inhomo- 
geneous frequency shift, i.e., w„(i) = f(i), which means 
that pi(i) = on the basis of Eq. (|10d[) . there exist two 
cases to be discussed: 

i) If p\H) 7^ in which we havepi(i) = pi = const 7^ 0, 
then this means that tp n it) is still a linear function of the 
discrete space n, i.e., (p n {t) = pin+po{t). In this case, the 
variable function r(t) and the inhomogencous frequency 
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FIG. 1. (color online). Profiles of the coefficients of the generalized ALH lattice with variable coefficients given by Eq. ([3]) vs 
time for the parameters given by Eq. (|17|) with parameters 70 = ci = ox = 1. (a) the coefficient vi(t) given by Eq. (|13a|) of the 
first degree term of the inhomogeneous frequency shift v„(t) given by Eq. (|12b|) . (b) nonlinearity g(t) given by Eq. (|12a|) . and 
(c) the gain or loss term j(t). 



shift v n (t) are given by 



A. Nonautonomous discrete one-rogon solution 



r(t) = 



Iq a{s)ds 



a cos 



i(pi) + fesin(pi) ' 

aa{t) p {t) 



acos(pi) + 6sin(pi) 



and p(t), g(t) are same as Eqs. and (|12a|) ; 

ii) If p±(t) = 0, i.e., tan-^&afi) - aj3(t)]/[aa(t) + 
bf3(t)]} = 0, which means that the relation for the coeffi- 
cients in Eq. ((3|), a(t) = (a/b)/3(t), is required and f n (t) 
is only a functions of time, i.e., ip n {t) = Po(t), then the 
variable function r(t) and the inhomogeneous frequency 
shift v n (t) are given by the form 



Tit) = - 

a 



a(s)ds 1 



i(t) =«(*)- 



Po(*) 



and p(t), g(t) are same as Eqs. (jllaj) and (|12aj) . where 
j(t), a(t) and Po(t) are free functions of time, and 
a, b, po are all free parameters. 



V. NONAUTONOMOUS DISCRETE ROGON 
SOLUTIONS AND INTERACTION 

In general, we have a large degree of freedom in choos- 
ing the coefficients of similarity transformation © and 
Eq. ([3]). As a consequence, we can obtain an infinitely 
large family of exact solutions of the generalized ALH lat- 
tice with variable coefficients given by Eq. @ in terms 
of exact solutions of the discrete Hirota equation §5§ and 
transformation (J8)). In particular, if we consider discrete 
rogon solutions of Eq. ([5]) as seeding solutions, then we 
can obtain many types of nonautonomous (including ar- 
bitrary time-dependent functions) discrete rogon (rogue 
wave) solutions of Eq. ([3]). As two representative exam- 
ples, we consider the lowest-order discrete rogon solutions 
of Eq. §5§ as two examples [33| to study the dynamics of 
nonautonomous discrete rogon solutions of Eq. J3]) . 



Firstly, based on the similarity transformation ([8j and 
one-rogon solution of the discrete Hirota equation ([5]), we 
present the nonautonomous discrete one-rogon solution 
(also known as the first-order rational solution) of Eq. ([3]) 
in the form 



*, ( i 1) (0 = PoVMCxp { - J j{s)ds + i[(p n (t) + M^} 



4(1 + /i) [l + Mpsja 2 + b 2 T{t)} 
1 + 4/m 2 + 16/x 2 (l + p)(a 2 + b 2 )r 2 (t) 

where the part phase <p n (r) is defined by 



(15) 



n (r) = 2r(i) {l+p)^/a 2 + b 2 -a -ntan -1 (6/a) , (16) 

p is a positive parameter, the variable r{t) is given by 
Eq. (jllcp . and the phase ip n (t) = Pi(t)n + po(t) with 
Pi(t) given by Eq. (jTTbJ) and po(t) being an arbitrary 
diffcrcntiable function of time. 

To illustrate the wave propagations of the obtained 
nonautonomous discrete one-rogon solution (|15[) . we can 
choose these free parameters in the form 

a(t) = c\ sin(2t), (3{t) = C2 cos(t), 

7(t) = 7 osin(t)cos 2 (0, (17) 
a = b = p = p = 1, 

where 70, are constants. 

Figure [T] depicts the profiles of the coefficient v-y{t) = 
[a(t)P(t)-a(t)/3(t)]/{2[a 2 {t)+f3 2 (t)]} of the first degree 
term of the inhomogeneous frequency shift v n (t) given 
by Eq. (|13a[) . the attractive intersite nonlinearity g(t) 
given by Eq. (|12a|) . and the gain or loss term j(t) vs 
time for the parameters given by Eq. (fl~7l). The evolu- 
tion of the intensity distribution for the one-rogon so- 
lution given by Eq. (jT5j) is illustrated in Fig. [5] for pa- 
rameters 70 = Ci = C2 = 1. The discrete rogue wave 
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FIG. 2. (color online). Profiles of nonautonomous discrete 
one-rogon solution (|15|) of the generalized ALH lattice with 
variable coefficients given by Eq. ([3]) for the parameters given 
by Eq. (fT7)l with 70 = Ci = C2 = 1. (a) the intensity 
distribution (t)\ 2 with max (njt) I*! 1 ' (t)| 2 4.7, (b) the 
density distribution ^^'(t)! 2 , (c) the intensity distributions 
l^n^t)! 2 for t = 0, 0.2, which means that the amplitude de- 
creases as time increases, and the peak falls the lower position 
after time exceeds about 1. 
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FIG. 3. (color online). Profiles of nonautonomous discrete 
one-rogon solution (|15|) of the generalized ALH lattice with 
variable coefficients given by Eq. ((3| for the parameters given 
by Eq. Q7J with 70 = 1, Ci = 0.01, and c 2 = 0.02. (a) 
the intensity distribution {^^(t)] 2 with max ( „ it) I*! 1 ' (f)] 2 S 
43.2, (b) the density distribution \^„\t)\ 2 , (c) the intensity 
distributions (t)| 2 for t = 0, 3, which means that the 
amplitude decreases as time increases, and the peak falls the 
lower position after time exceeds about 70. 



solution is localized both in space and in time, thus re- 
vealing the usual discrete 'rogue wave' features. However 
if we fix the coefficient 70 = 1 of the gain or loss term and 
adjust the coefficients c\ = 0.01, C2 = 0.02 of the tun- 
nel couplings a(t) and j3(t) given by Eq. (| 1 T[) . then the 
evolution of the intensity distribution for the one-rogon 
solution is changed (see Fig- [3]) , and it follows from Fig. [3] 
that the discrete one-rogon solution in this case is local- 
ized in space and keep the localization longer in time 
than usual rogue waves (see HH). Moreover, it follows 
from Figs. [2{c) and [3jc) that the amplitude of the dis- 
crete one-rogon solution decreases as time increases, and 
the amplitude in Fig. HJc) is decreased faster than one in 
Fig. |3I c) as time increases. 



B. The interaction between nonautonomous 
discrete rogon solutions 

Here we consider the interaction between nonau- 
tonomous discrete rogon solutions. To do so, we ap- 
ply a second-order rational solution of the discrete Hi- 
rota equation ((5]) to the similarity transformation ([8]) 
such that we can obtain the nonautonomous discrete two- 
rogon solution of Eq. © in the form 

^nKt) = poV^CXp j( S )ds + i[tp n {t)+ip n (T)}\ 

12(l + / ,)K 2) (r) + Jv /^Qi 2) (r)]" 
X 1 (2), . 18 

which displays the interaction between nonautonomous 
discrete rogon solutions, where [i is a positive parameter, 

the functions Vu\t), Q« (t) and Hn\r) arc all poly- 
nomials of space and time given by 

vL 2) (t) = 5^ 2 + 6(JK + 2^ + 3)^ + 
+(6-4^)^-3(4^ + 1), 

Q^\t) = ^ 2 + 2(^+l)^ + ^+ 2 

-(16/i + 6)^-3(8^ + 5), 

H { n ] (r ) = ,r 3 + 3 {JV + 8/x + 9) .J 2 + 3 („+ 2 - 6^+ 
-16^+48^ 2 + 72/i+33)^+.yK 3 
+(3 - 8/i)^K 2 + (27 + 24/i + 16/^)^+9, 

where we have introduced jV = Afin 2 and 8? = 16/z 2 (l + 
^)(a 2 + 6 2 )r 2 , t = r(t) is given by Eq. (jllcj) . the part 
phase <p n (t) = Pi(t)n+po(t) with pi(i) given by Eq. ()llbj) 
and po(t) being an arbitrary differentiable function of 
time, and the part phase (p n (j) is given by Eq. (|16l) . 

Similarly, we can choose these free parameters given 
by Eq. (|17|) for the nonautonomous discrete two-rogon 
solution given by Eq. (fT8| except for fi = 1/16. Fig- 
ures @] and [5] depict the evolution of intensity distribu- 
tion for the interaction between nonautonomous discrete 
rogon solutions (discrete two-rogon solution) given by 
Eq. (|18|) for different parameters 70, c\, and C2, respec- 
tively. Moreover, we find that the interaction between 
nonautonomous discrete rogue waves is elastic. Figure 0] 
shows that the nonautonomous discrete two-rogon solu- 
tion is localized both in space and in time, thus revealing 
the usual discrete 'rogue wave' features for the chosen 
parameters 70 = 1, c\ = 2, C2 = 1, but if we fix the co- 
efficient 70 = 1 of the gain or loss term and adjust the 
coefficients c\ = 0.2, C2 = 0.1 of the tunnel couplings a(t) 
and P(t) given by Eq. (JT7J) , then the evolution of the in- 
tensity distribution for the discrete two-rogon solution is 
changed (see Fig. [S]), and it follows from Fig. [5] that the 
nonautonomous discrete two-rogon solution is localized 
in space and keep the localization longer in time than 



6 



usual rogue waves. Moreover, it follows from Figs. |4jc) 
andlSfc) that the amplitude of the discrete two-rogon so- 
lution decreases as time increases, and the amplitude in 
Fig. Utc) decreases faster than one in Fig. [SJc) as time 
increases. 




FIG. 4. (color online). Profiles of the interaction between 
nonautonomous discrete rogon solution given by Eq. {THJ of 
the generalized ALH lattice with variable coefficients given 
by Eq. @ for the parameters given by Eq. (|17p with 70 = 
1, ci = 2, and C2 = 1. (a) the intensity distribution 
|*l 2) (£)| 2 with max (n , t) |*L 2) (f)| 2 1.65, (b) the density dis- 
tribution |^l 2) (t)| 2 , (c) the intensity distributions |^l 2) (t)| 2 
for t = 0, 0.6, which means that the amplitude decreases as 
time increases, and the peak falls the lower position after time 
exceeds about 3. 
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FIG. 5. (color online). Profiles of the interaction between 
nonautonomous discrete rogon solution given by Eq. 1)181) of 
the generalized ALH lattice with variable coefficients given 
by Eq. © for the parameters given by Eq. (|17p with 70 = 
1, ci = 0.2, and C2 = 0.1. (a) the intensity distribution 
|*l 2) (£)| 2 with max (n , t) |*L 2) (£)| 2 2.25, (b) the density dis- 
tribution |^l 2) (t)| 2 , (c) the intensity distributions |^l 2) (t)| 2 
for t = 0, 1.8, which means that the amplitude decreases as 
time increases, and the peak falls the lower position after time 
exceeds about 65. 



VI. CONCLUSIONS 

In conclusion, we have studied nonautonomous dis- 
crete rogon solutions and their interaction in the general- 
ized Ablowitz-Ladik-Hirota lattice with the varying tun- 
nel coupling, intersite nonlinearity, inhomogeneous fre- 
quency shift, and gain or loss term given by Eq. ([3]) on 
the basis of the differential-difference similarity reduction 
J8]). Wc found its some nonautonomous discrete rogon 
solutions when the intersite nonlinearity g(t) (the inho- 
mogeneous frequency shift v n (t)) is related to the gain or 
loss term j(t) (the tunnel coupling A(t) = a(t) + i(3(t)) 
(see Eqs. (|12a[) and (|12b|) ). This denotes that only two 
coefficients (i.e., ^(t) and A(t) = a(t) + if3(t)) are left 
free. 

In particular, we have studied wave propagations of 
nonautonomous discrete rogon solutions and interaction 
for some chosen parameters, which exhibit complicated 
rogue wave construes. For the given periodic gain or 
loss term, when the amplitude of the tunnel coupling 
coefficient between sites decreases, these nonautonomous 
discrete rogon solutions are localized in space and keep 
the localization longer in time, which differ from the usual 
discrete rogue waves of nonlinear discrete equations (e.g., 
the AL lattice and the discrete Hirota equation) [33[ . 

Moreover, nonautonomous discrete rogon solutions 
and interaction may provide more documents to fur- 
ther understand the physical mechanism of discrete rogue 
wave phenomena. The approach may also be extended to 
other discrete nonlinear lattices with variable coefficients 
for studying their discrete rogue wave solutions and wave 
propagations. 
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